Image deformation in field ion microscopy of faceted crystals 
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We perform detailed numerical simulations of field ion microscopy images of faceted crystals and 
compare them with experimental observations. In contrast to the case of crystals with a smooth 
surface, for a faceted topography we find extreme deformations of the ion image. Local magnification 
is highly inhomogeneous and may vary by an order of magnitude: from 0.64 to 6.7. Moreover, the 
anisotropy of the magnification at a point located on the facet edge may reach a factor of 10. 
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I. INTRODUCTION 



In recent years field ion microscopy (FIM)^'- has 
found a number of applications connected with ther- 
mally faceted crystals^*^. The most prominent exam- 
ple is the research on the fabrication of atomically 
sharp electrodes^iii^iSiiiSiii and their use as practical 
electror*i^ and ion^- point sources, including laser-driven 
femtosecond electron source a^'^d^ . In these studies FIM 
is often used to determine the electrode shape and the 
degree of its sharpness. Another example of the applica- 
tion of FIM to faceted crystals is its use in surface science 
studies of the faceting process itself, in its steplike and 
global forroi^iiiii^, and in the study of the equilibrium 
crystal shape^^. Worth mentioning are also FIM studies 
on single atom diffusion^*^ - although not carried out on 
thermally-faceted surfaces, they also encounter the prob- 
lem of large image distortion at the facet edge. 

It is well established that image magnification in FIM 
is not uniform-. This constitutes a problem in the inter- 
pretation of the FIM images, and in the past researchers 
have tried to devise methods which would allow to take 
into account these effects. Both analytical and numerical 
calculations of the electric field and ion trajectories have 
been carried out^'^^'^^. However, none of these calcula- 
tions apply to faceted crystals, which are recently most 
often examined in FIM. Our work is aimed at filling this 
gap. Faceted crystals, with their atomically sharp edges, 
produce very large local variations of the intensity and 
direction of the electric field. This produces huge deflec- 
tions in ion trajectories, and, consequently, huge distor- 
tions in the resulting microscopic image. 

In this paper we perform detailed numerical calcula- 
tions of the ion trajectories in the FIM of faceted crystals. 
This leads to surprising insights into the interpretation 
of field ion images of such crystals. Where possible, we 
confront our numerical results with available experimen- 
tal data concerning image distortion. 



II. METHOD OF CALCULATION 

In this section we discuss the numerical approach to 
simulation of field ion microscope images of faceted sur- 
faces. First, we find the electric field distribution in a 
model of FIM by solving Laplace's equation, and then 
the ion trajectories are obtained form classical equations 
of motion. 



A. Construction of the model of a faceted crystal 

We will use specimens generated earlier in Monte Carlo 
simulation of adsorbate-induced faceting on curved sur- 
faces of bcc metalS^. A typical specimen has spherical 
shape with faceted area around the the [111] pole limited 
by the inclination angle Q — 20° . The faceted region of a 
specimen has a form of a single pyramid with {112} faces 
or it contains step-like {112} mircofacets. We also study 
truncated pyramids and a pyramid with double edges. 

We are especially interested in answering the question 
how the morphology of faceted surface affects the elec- 
tric field distribution near the emitter surface. Moreover 
we want to examine the local magnification on images of 
faceted surfaces. Thus the FIM model should take into 
account the atomic roughness of the tip surface and as- 
sure long enough distance between the emitter and the 
screen, or in other words, the model should span many 
length scales. Such conditions cause big difficulties in nu- 
merical solution of Laplace equation on a discrete mesh. 
On the one hand, the size of the mesh should be equal 
to the distance between the tip and the screen ( ~10 cm 
in real experiment). On the other hand, the mesh should 
have resolution smaller than 1 A near the tip. In this pa- 
per we demonstrate that solution of Laplace's equation, 
in a case when the distance between the emitter and the 
screen is million times larger than the mesh resolution, 
is possible assuming the following approximations: 

• The spherical emitter with faceted region. Atoms 
in the faceted region are represented by the trun- 
cated octahedrons - the Wigner-Seitz cells of the 
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FIG. 1: 3D view of the apex of a pyramid formed of the 
truncated octahedrons. 



bcc crystal structure (see Fig. [T]). The remaining 
part of the emitter is approximated by a smooth 
sphere of radius Rq. 

• The spherical screen of radius i?i which is about 
three hundred times grater than Rq. 

Let us comment the the choice of truncated octahe- 
drons to represent atoms in the FIM model. We con- 
sidered 3 possibilities: (i) point like atoms and the sur- 
face represented by triangles, (ii) hard sphere model, 
(iii) each atom represented by the truncated octahedron 
- the Wigner-Seitz cell of a bcc crystal structure. All 
these cases were preliminarily tested in simulations and 
we have chosen the representation by truncated octahe- 
drons. The simplest representation by surface triangu- 
lation does not accurately reflect the atomic roughness 
of (111), (211) and (110) faces. The hard sphere model 
seems to be the best one, but it is difficult to use. If one 
takes the radius from the close-packed bcc structure, the 
surface contains plenty of holes between the spheres. On 
the other hand, choice of larger sphere radius leads to 
the reduction of the atomic roughness. Representation 
of atoms by truncated octahedrons ensures that the bcc 
crystal faces have proper symmetry and atomic rough- 
ness, and it is easy to use in numerical calculation. 



B. Calculation of the electric field 

To calculate the electric field distribution in the FIM 
model, first one has to solve Laplace's equation (LE) in 
the space between the emitter and the screen^^'^^. 



AV{r) = 0, 



(1) 



where V{r) is the electric potential at the point r. It is 
assumed that V = Vq a.t all points of the emiter surface 
and the V — Vi at all points of the screen. Then the 
electric field is obtained: 



E(r) = -VV{r). 



(2) 



To solve LE, Eq. in the space between the emit- 
ter and the screen, finite element method (FEM) is ap- 
plied, because it is well-suited for simulation of geometri- 
cally complicated domains^'*'^^. All FEM related calcula- 
tions are performed using GetfemH — h packageSS, whereas 



tetrahedral meshes are generated using TetGen library^!. 
In our application of the finite element method the do- 
main is divided into a number of tetrahedrons and 
Lagrange-type interpolation functions are used as the ba- 
sic functions. For such mesh the electric field is constant 
in a tetrahedron. Hence, to assure appropriate accuracy 
of E, especially near the emitter surface, the linear size 
of a tetrahedron should be much smaller than the lattice 
constant a. Direct computation does not allow for such 
mesh refinement, so we overcome this difficulty in the 
following way. First, we combine the numerical solution 
of LE with the analytical one. It follows from our pre- 
liminary calculation that at distances r from the emitter 
greater than 5Ro the electrostatic potential V{r) for the 
faceted emitter is practically the same as for the spherical 
emitter. 

It is easy to show that the electrostatic potential Vg 
for the spherical emitter is given by the formula 



Vsir) 



ViRi-VoRo , {Vo-V,)RoRil 



Rl — Rq Ri — Rq 

and applying Eq. ^ gives the electric field Eg 



E5(r) = 



{Vq - Vi)RoRi r 
Rl - Rq r3 



(3) 



(4) 



A numerical solution of LE is limited to a subdomain 
n' of spherical shape with a radius i?num- We apply 
Dirichlet boundary conditions to surfaces in $7' with V — 
Vq on internal (emitter) surface and V = Vs(i?nuni) on 
the external surface. In the space outside the fl' , i.e., 
for distance r > i?num, the analytical solution given by 
Eq. ^ and Eq. (g]) will be used. 

To obtain a numerical solution of LE in the subdomain 
O' with appropriate accuracy for FIM images simulation, 
we use 3-step procedure described below: 

1. The domain f2' is divided into n spherical layers - 
subdomains Qi, i — 1, 2, . . . , n, each with its own 
mesh of tetrahedrons. Then the LE is solved in the 
whole rt' to get values of V on boundaries of each 
subdomain. 

2. In each subdomain Qi the LE is solved again on 
a finer mesh of tetrahedrons applying Dirichlet 
boundary conditions to fli surfaces where the elec- 
tric potential is known from solution obtained in 
the previous step. 

3. To calculate the electric field at a point r within 

(a) A local subdomain tor is constructed by col- 
lecting all tetrahedrons with distance from r 
smaller than i?ioc. 

(b) A new mesh is created by filling with newly 
generated tetrahedrons the space inside the 
boundary. There is a special, regular tetrahe- 
dron Tr with center located at r 
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(c) Using this mesh, LE is solved imposing V val- 
ues from the ujj^ boundary as Dirichlct condi- 
tion. 

(d) Finally, the electric field is obtained as minus 
gradient of V calculated on Tr 

Typical technical data used in our calculation: the ra- 
dius of emitter Rq — 268 A, assuming the lattice con- 
stant for tungsten: a=3.16 A, the radius of 51' domain 
^num — 33i?o, the radius of the screen Ri = 267 Rq, the 
number of subdomains n — 9, the number of tetrahedrons 
used in whole W for the 1st step is Nt = 4728000 and 
for the 2nd step in each subdomain O^, Nt > 2.6 * 10^. 
The size of the special tetrahedron Tr is ^ in the vicinity 
of emitter surface. Hence, the ratio of Ri to size of the 
special tetrahedron is over 10^. 

It is worth to emphasize an important role of the third 
step in our procedure of calculating the electric field by 
using a local subdomain ujr- Omitting this step leads to 
an incorrect solution even inQi subdomain, although Qi 
has the smallest linear size 1.07i?o) and large number 
of tetrahedrons in its mesh Nt = 3.4 x 10^. 

We performed several tests to verify presented method 
of solving LE. The most important one concerned cal- 
culation of the electric field for an ideal spherical emit- 
ter because the analytical solution is known in this 
case - Eq. ([4]). The obtained numerical results are in 
good agreement with analytical solution. In the exam- 
ined distance range up to r = 8i?0j the relative error 
\E{r) — Es{r)\/Es{r) is smaller than 5*10~^ and the de- 
viation of E direction from the radial direction is smaller 
then 0.1°. We also obtained positive results for check- 
ing the symmetry of the electric field distribution on the 
different emitter shapes: pyramidal tip and bcc spherical 




FIG. 2: An example of ion trajectories (continuous lines) and 
electric field lines (dashed lines) starting from the edge of a 
pyramid. In this image the length of these lines is limited to 
0.5Ro. 



voltage, mass and charge of the ion. These quantities 
affect unit of time tQ. To solve Eq. ^ we apply the 
velocity Verlet algorithm - commonly used in molecular 
dynamics simulations^S., with a time step O.lio- 

It is obvious that the direction of ion trajectory is not 
coincident with the electric field lines (see Fig. [2]) even in 
the case of the initial ion velocity equal to zero. Due to 
the curvature of the electric field lines, instantaneous ion 
velocity at r is deflected from E(r) up to large distances 
from the tip surface. 



D. Calculation of a local magnification 



C. Calculation of the ion trajectories 

Motion of an ion of mass M and charge q under the in- 
fiuencc of force F(r) — gE(r) is described by the equation 
of motion 



M 



qE(r). 



(5) 



It is assumed that the ion is formed with zero initial 
velocity at a small distance (a-^) from a surface atom. 
It is convenient to express quantities occurring in Eq. ([5]) 
in reduced units. Choosing lattice constant a as unit of 

length, (Vo - Vi)/a as unit of E and to — \J q(Vo-v ' i^ 
as unit of time, the equation of motion Eq. ^ can be 
expressed in the form 



(6) 



where tilde denotes quantity in reduced unit. Thus, the 
trajectories in reduced units do not depend on the applied 



Knowing the full ion trajectories we can investigate 
magnification not only on the distant microscope screen 
but also on virtual screens at small distances from the 
emitter. Analysis of FIM images on virtual screens is 
used to determine the minimal emitter - screen distance. 
Magnification / is defined here as 



d' 



(7) 



where d is the distance of a pair of points on the specimen 
and d' is the corresponding distance on the screen. 

Another important quantity is the local magnification 
// of distance d defined as 



d' 



fi 



fav 



(8) 



where fav stands for average magnification. The average 
magnification can be replaced by magnification of a large 
distance, eg. the distance between ends of pyramid edges. 

We found that local magnification of faceted crystal 
has a long-range dependence on the distance r from the 
tip. Using the tip with curvature radius Rq « 27 nm 
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we obtained that the local magnification reaches 73% 
of its final value /f at r = Rq. At larger distances 
M5Ro) - 0.9/f, /z(10i?o) « 0.94/f, and //(30i?o) « 
0.98/f. Hence, to control the error of /;, the screen - 
tip distance in the numerical model should be at least 
30 times greater than the tip curvature radius, which is 
much greater than the distances used in previous numer- 
ical studies (see egj ^^'^^ ). 

III. RESULTS AND DISCUSSION 

A. A perfect, single atom pyramid 

We begin our calculation with the simple case of a 
perfect 3-sided pyramid pointing in the [111] direction, 
formed by three densely packed crystal facets: (211), 
(121), (112). The crystallographic direction of the pyra- 
mid edges is (311). We assume that all atoms are in the 
BCC lattice positions and that the pyramid is ended by a 
single atom. This implies that that the number of atoms 
in the succcesive (111) planes is 1,6,15,28.... The as- 
sumed sample configuration is shown in Fig.[3ja) and (b). 
After calculating the resulting electric field E{r) between 
the sample and the microscope screen, we map the sur- 
face atoms onto the microscope screen by following the 
classical charged-particle trajectories (ions or electrons) 
- the solution of Eq. [S] The result is shown in Fig. OJc). 
The calculated image shows some striking features. 

1. Of all the straight lines seen on the sample 
[Fig-E^b)], only the pyramid edges remain straight 
lines on the image [Fig. [3fc)]. 

2. The magnification of the image is highly inhomoge- 
nous. There is a great enhancement of the mag- 
nification at the apex of the pyramid. The local 
magnification factor at the apex is 4.5. 

3. Along the pyramid edge, the local parallel mag- 
nification varies from 2.8 near the apex, reaching 
a minimum (0.64) near the center of the edge (de- 
magnification), increasing to a value of 1.5 near the 
pyramid base. 

4. Along the pyramid edge, the local transverse mag- 
nification is 7.4 near the apex, quickly decreasing 
to a constant value of 6.7. 

5. As seen from the comparison of features [3] and 
Ul, the magnification of the image is exceedingly 
anisotropic - the image is extremely stretched in 
one direction. In the middle of the pyramid edge, 
the anisotropy is 6.7/0.64 ~ 10. 

It would be very interesting to confront the above re- 
sults with experiment. However, there is no easy way 
to calibrate an FIM image in terms of real distances. It 
is possible with careful FIM of samples with low aver- 
age radii, where individual atoms are resolved, but as 




FIG. 3: A single atom pyramid, (a) A three-dimensional 
visualization of the sample, (b) A projection of the surface 
atoms of the sample onto a plane, (c) Calculated projection 
onto the microscope screen following the charged-particle tra- 
jectories in the electric field. 



the average radius increases, the microscope resolution 
decreases and one looses the natural reference length of 
the lattice constant. However, observation [51 appears to 
explain a peculiar feature of FIM images of pyramids - 
the images of atoms at the edges are often not circular, 
but stretched in the direction perpendicular to the edges 
- see Fig. |31 

The image of atoms shown in Fig. [3fc) does not corre- 
spond directly to the experimentally observed FIM im- 
age, because in the real microscope only the most pro- 
truding atoms are visible. One of the reasons is that the 
protruding atoms, which can be thought of as "sharp" 



FIG. 4: A FIM image of oxygen-covered tungsten pyramidi^. 




FIG. 5; Atoms which satisfy the ionization criterion (above a 
treshold field): (a) the sample, (b) the simulated FIM image. 



features on the otherwise sniooth surface, cause local en- 
hancement in the electric field. Imaging requires field ion- 
ization, which takes place for electric field magnitudes ex- 
ceeding the ionization treshold. Experiments show that 
for the crystal shape considered here, only the edges of 
the pyramid are seen in the FIM image. For this reason, 
for the realistic simulation of the image, we select the 
atoms where the electric field above the atom exceeds a 
certain treshold [Fig.tSJa)]. In the last step we calculated 
the ion trajectories starting from the seletected atoms 
and extending up to the microscope screen to generate a 
simulated FIM image [Fig. [SUb)]. 



FIG. 6: Visualization of a truncated pyramid (7 atomic (111) 
layers removed). 



B. Truncated pyramids 

The exact shape of the {211}-faceted crystal can be 
controlled by careful temperature treatment. It is pos- 
sible to obtain pyramids with various degrees of vertex 
truncation or vertex rounding, which can be quantified 
by measuring the distance d between the pyramid edge 
end pointsi^. However, such measurement in FIM suffers 
from errors caused by the inhomogenity of the local mag- 
nification. If one assumes constant magnification in FIM, 
the measurements of d are always overestimated - that is, 
the pyramids are actually more "sharp" than they appear 
to be in the microscopic image. To calculate the correc- 
tion for the measurements carried out in^, we calculated 
a series of images of truncated pyramids. The height of 
the single atom pyramid was 14 geometrical atomic lay- 
ers; now we remove 1,2,3 etc, obtaining an atomic con- 
figuration shown in Fig. [5] Three examples of resulting 
images are shown in Fig. [71 A striking feature is the 
shape deformation of the truncated vertex area, which 
has a triangular boundary on the assumed sample, but 
almost a circular boundary on the simulated image. In 
FIM experiments, the truncated area is not well resolved, 
but the boundary of the truncated vertex indeed often 
appears curved outword, as shown in Fig. [51 As could be 
expected from the results presented in Sect. lIIlSl the lo- 
cal magnification in the central part of the image is high. 
For one atomic layer removed, the area of the truncated 
vertex increases 8.0 times; for 2 layers - 4.8 times; for 
5 layers - 2.1 times. The corresponding linear magni- 
fications, proportional to the square roots, are 2.8, 2.2 
and 1.4, respectively. As the total height of the pyramid 
was 14 layers, the local magnifications calculated above 
correspond to height truncation of 7%, 14% and 36%. 

Applying similar considerations as above to the mag- 
nification of the distance d between the edge end points, 
it is possible to correct previous work on the dependence 
of the equilibrium crystal shape on the temperatureiS,, as 
shown in Fig. [51 Note that without the correction based 
on this work, the amount of truncation of the pyramid is 
highly overestimated. 

In Fig. [ini we demonstrate the simulated image of the 
(111) surface formed by truncating half of the height of 



FIG. 7: Truncated pyramids with (a) 1, (b) 2, and (c) 5 atomic layers removed. The corresponding simulated ion images are 
shown in (a')-(c') 
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FIG. 8; A FIM image of a truncated pyramid. The arrows 
point to the edge end points. 



the {211} pyramid. In the experiment, at this pyramid 
size the indywidual atoms would not be resolved; how- 
ever, in principle all the (111) lattice sites could be visu- 
alized using single atom adsorption on the (111) plane, 
similarly as in the work of Antczak and Ehrlich^°. In 
this way they have observed huge deformations of the 
image of an atomically flat crystal facet: curving of the 
atomic lines and an increase of the magnification near the 
boundary of the facet (Fig. 2 in Refi^). Both of these 
features are present in the result of our simulation shown 
in Fig. I10[ further confirming the validity of our model. 
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FIG. 9: The corrected dependence of the equilibrium crystal 
shape on the temperature for oxygen-covered tungsten^^. 



C. A steplike-faceted crystal (hill-and-valley 
faceting) 

In the experiment, occurence of global faceting, where 
the shape of the faceted crystal is convex, such as de- 
scribed in Sect. IIII Al and IIII Bl is only possible under 




FIG. 10: (a) The atomic arrangement on the (111) plane of 
the truncated pyramid in Fig. [6] (c) Calculated projection 
onto the microscope screen following the charged-particle tra- 
jectories in the electric field. 



special conditions: small crystal, high annealing tem- 
perature, and high desorption temperature of the adsor- 
bate. If these conditions are not fulfilled, one obtains a 
steplike-faceted crystalS'liS. This is due to the kinetic lim- 
itations on the surface diffusion of the crystal material. 
For this reason we consider here a steplike-faceted sam- 
ple, as shown in Fig. 1111 The assumed atomic configura- 
tion closely corresponds to surface topography observed 
in the experiments (see e.g. RefJ^). 

Comparison of the simulated and real FIM image 
[Fig. fllTc) and (d)] shows an overall agreement, but also 
a discrepancy in one aspect. In the experiment, the six 
(311) step edges are not pairwise parallel. One possible 
reason is that our model sample has a different topogra- 
phy than the real sample near the boundary of the faceted 
region. Another possibility is that the real atomic con- 
figuration is not as perfect as assumed in the calculation 
[Fig. [TlTa)]. Note the high variation of the local mag- 
nification along the step edge in the simulated image. 
Unfortunately, this effect cannot be verified in the real 
image [Fig. [TTT d)] due to the lack of single-atom resolu- 
tion. 

Another atomic configuration, which is interesting for 
comparison with experiment, is shown in Fig. 1121 This is 
significant for the case of palladium-induced faceting of 
tungsteni^iS^. In these experiments, palladium forms a 
pseudomorphic physical monolayer on the sample. The 
(311) edges of pyramids/steps are truncated, with one 
atomic row removed, as shown in Fig. 1121( a). Compari- 
son of the actual atomic arrangement [FigfT2yb)] with the 
simulated FIM image at the microscope screen [FigfT^c)] 
reveals a serious problem in the interpretation of the FIM 
image in such a case: it is very difficult to infer from the 
image on the screen [FigfT2lc)] the actual configuration of 
the atoms [FigfTWb)]. The image observed on the screen 
is very misleading, suggesting the presence of three se- 
ries of four almost equally-spaced parallel edges. This 
constituted a serious problem in the interpretation of the 
ion images of Pd/W faceting. Our present calculations 
finally fully confirm the interpretation of the ion images 
assumed in previous work^^'— , illustrated here in Fig. 1131 




FIG. 11: A steplike-faceted crystal, (a) 3D visualization of 
the assumed atomic configuration. (b),(c) The positions of 
the image atoms on the sample and on the screen, (d) Com- 
parison with a real image. 



Calculations yield a very large magnification anisotropy: 
in the middle of the facet edges, the local transverse mag- 
nification is 3.1, while the local parallel magnification is 
0.74. This results in the image magnification anisotropy 
of 4.2 at this point - the FIM image is stretched over 
4 times in the direction transverse to the facet edges. 
In the experiment, the local transverse magnification is 
probably even larger, a rough estimate yields a factor 
of 7-8 (compare Fig. [T3|) . However, this number is only 
an estimate, as there is no direct way to measure local 
magnification in the FIM experiment. 
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FIG. 12: (a),(b) The assumed atomic configuration and (c) 
the simulated ion image, corresponding to Pd/W faceting. 



IV. CONCLUSION 

We performed detailed trajectory calculations in the 
FIM of {211}-faceted crystals (the case often encountered 
in research), and compared them with available experi- 
mental data. We found that the experimentally obtained 
FIM images of faceted crystals are often very misleading, 
and one has to rely on numerical calculations to fully 
account for huge local deformations of the microscopic 
image. Our main conclusions can be summarised as fol- 
lows: 

1. In the consideration of faceted crystals, it is neces- 
sary for the model of FIM to take into account the 
atomic roughness of the tip surface and, at the same 
time, the large distance from the tip to the screen 




FIG. 13: (a) FIM image of a faceted palladium-covered tung- 
sten surface^^. (b) The position (311) edges in image (a), 
(c) Approximate actual position of the edges on the crystal 
surface. The distances are estimated from FIM and STM 
experiments^^. 

- the model must span many length scales. Our so- 
lution of Laplace's equation demonstrates that the 
distance between the tip and the screen in the FIM 
model should be at least 30 times greater than the 
tip curvature radius in order to calculate the local 
magnification correctly. 

2. A combination of the finite element method with an 
analytical solution allows for accurate calculation of 
the electric field in the space between the tip and 
the screen. 

3. The FIM images of faceted crystals exhibit very 
large variations of the local magnification. We have 
found local magnification factors as high as 6.7 and 
as low as 0.64 (local de-magnification). 
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4. As one moves along the pyramid or step edge, the 
paraUel magnification varies by a factor of 2-4, 
reaching a minimum near the middle of an edge. 

5. The images of pyramid edges or step edges are ex- 
ceedingly stretched in the direction perpendicular 
to the edge. The anisotropy of the local magnifica- 
tion, defined as the ratio of local magnification in 
two ortogonal directions, reaches a factor of 10 (for 
a perfect {211} pyramid). 

6. As a consequence of [5], an image of a truncated 
facet edge mimics the image of two distant parallel 
edges. Great caution is required in the interpreta- 
tion of this kind of FIM images. 

7. The image of a single flat facet is distorted in such 



a way that the local magnification increases as one 
moves toward the facet edge. The atomic rows do 
not appear straight on the image. Although we 
have calculated only the case of a (111) facet on top 
of a {211}-pyramid, comparison with experiment 
suggests that this conclusion is generally true for 
isolated flat facets. 



8. Truncated {211} pyramids appear significantly 
more truncated in the FIM image than they are 
in reality; this effect is most pronounced for small 
truncations. Numerical calculations are necessary 
to estimate the size of a {211}-faceted tip. 
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